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Abstract. We show that fractional exclusion statistics is manifested in general in in- 
teracting systems and we discuss the conjecture recently introduced (J. Phys. A: Math. 
Theor. 40, F1013, 2007), according to which if in a thermodynamic system the mutual 
exclusion statistics parameters are not zero, then they have to be proportional to the di- 
mension of the Hilbert space on which they act. By using simple, intuitive arguments, but 
also concrete calculations in interacting systems models, we show that this conjecture is not 
some abstract consequence of unphysical modelling, but is a natural-and for a long time 
overlooked-property of fractional exclusion statistics. We show also that the fractional ex- 
clusion statistics is the consequence of interaction between the particles of the system and it 
is due to the change from the description of the system in terms of free-particle energies, to 
the description in terms of the quasi-particle energies. From this result, the thermodynamic 
equivalence of systems of the same, constant density of states, but any exclusion statistics 
follows immediately. 
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1 INTRODUCTION 

Fractional exclusion statistics (FES) is a generalization of Bose and Fermi 
statistics, introduced by Haldane in Ref. [1], and with its thermodynamic 
properties calculated mainly by Isakov [2] and Wu [H] . The concept have been 
applied to a large number of systems of interacting particles [H El El El El El 
[T0| [TTl [T2l [T3l 111] and describes quasiparticles that exist in finite dimensional 
Hilbert spaces-in general, quasiparticles in a finite region of condensed matter; 
the Hilbert spaces are extensive, increasing proportionally to the size of the 
condensed matter region [T]. Quasiparticles that exist in these Hilbert spaces 
are called in general species-each subspace contains one species. Let us denote 
the Hilbert spaces by Tii, with i = 0,1, . . ., each of them containing Ni particles 
and Gi available states |15J. Then, by increasing the number of particles of 
one species, say Ni increases to Ni + 6Ni, the number of available states in 
any of the Hilbert spaces changes by 6Gj = —ajiSNi. The proportionality 
factors, aij, are called direct (when i = j) and mutual (when i / j) exclusion 
statistics parameters. 




Figure 1: In (a), TLi and TL2 are two Hilbert spaces defined by the ranges of 
quasiparticle energies (the energy levels in 7ii are not relevant, and therefore 
are not represented). In (b), extra particles are inserted into Tii, which changes 
the density of states in "^2 and therefore the number of allowed single-particle 
states. This is a manifestation of fractional exclusion statistics. 

The simplest examples are the Bose and Fermi statistics. For ideal 
bosons, aij = for any i and j, whereas for ideal fermions an = 1 for any i 
and aij = for any i and if i 7^ i. 

2 PROPERTIES OF THE MUTUAL EXCLUSION 
STATISTICS PARAMETERS 

The FES (other than Bose and Fermi statistics) is the result of interaction 
between quasiparticles. The finite dimensional Hilbert spaces involved, Tii, 
are formed of quasiparticle wavefunctions corresponding to eigenvalues con- 
tained in finite regions of the phase space. For example Isakov p] and Iguchi 
and Sutherland [TBI E] define the Hilbert spaces by the wavevector eigenval- 
ues, k, whereas Murthy and Shankar [10], Sen and Bhaduri [llj . Hansson, 
Leinaas, and Viefers [12] define the Hilbert spaces by the quasiparticle energy 
eigenvalues, e. In all the cases, the eigenvalues of the wavefunctions contained 
in any of the Hilbert spaces belong to a finite range (or a multi-dimensional 
volume); therefore the quasiparticle-quasiparticle interaction, which changes 
the density of states (DOS) in the phase-space, changes also the number of 
quasiparticle states in the range, leading to FES [18] (see figured]). 
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Let's take for example a system where the quasiparticle energies may be 
written asfinillllllllElliaiSnillH] 



1 

h = ei + '^VijUj + -ViiTii, (1) 
i=o ^ 

where are the energies of the noninteracting particles; The total energy of 
the system is -E = X^i^o "-j^j- We assume that the system in large enough-and 
the energy levels dense enough-to introduce the (quasi) continuous density of 
states, a{e) and write equation ([T]) as 

~e = e+ rV{e,e')n{e')a{e')de', (2) 
Jo 

where we replaced the subscripts i and j of Vij by the corresponding energy 
levels, e and e', respectively. In the variable e, the density of states is different 
from cr(e) and we shall denote it by a{e). To calculate a{e), we take a small 
interval, 5e, containing a{e) ■ 5e states, and transform it into the interval de, 
which, obviously, will contain the same number of states. Dividing the number 
of states by the energy interval and using equation ([2]), we find 

'7[e(e)] = mrr-^^ — (3) 

1 + /o ^^'yie'MeVe' + V{e, e' / e)a{e)n{e) 

In figure [2] we show an example. In some arbitrary, dimensionless units, with 
the population of the single particle energy levels shown in the left plot, the 
DOS (T(e) shown as the solid line in the middle plot, and constant interaction 
potential, V{e,e') = V, we calculate 5'[e(e)] (point line in the middle plot) 
and e (point line in the right plot)-for the concrete calculations of figure Owe 
choose V = I, a{e) = y^, and n(e) = [exp(e + 1) — 1]^^- 

To show how FES is manifested in the system, we split the quasiparticle 
energy axis into small intervals, [eo, ei ],..., [ej_i, e,], .. . (as shown in figure Ej), 
so that each interval contains large enough numbers of particles and available 
single particle states; we denote by A^(ej,ej+i) = Ni the number of particles 
in the interval [ei,ej+i], 

Ni = Nieu ei+i) = a{e)n{e) de = cT{e')n{e') de' , 

and, assuming that the energy intervals are small enough so that we can replace 
in the integral of equation ([2]) V{eM,^) by V{eM,Gi-i) for any e G (ej,ej_i], 
we write 

i-i 

= <^j + T. ^(^J' ei+i)- (4) 

i=0 
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Figure 2: The DOS on the quasiparticle energy axis, cj[(e)]-the point line in the 
middle plot-, and the quasiparticle energy, e-the point line in the right plot- 
for a bosonic single particle energy levels population, n(e) = [exp(e + l) — 
a DOS (T(e) = -y/e and a constant interaction potential, V{e,e') = 1. In the 
right plot we show both, e(e) (point line) and e(e) (solid line) for comparison. 

If we insert 5Ni particles into the interval [ej,ej+i] and we keep all the free 
particle energy levels ej(>i) unchanged, then we change the quasiparticle energy 
levels ej(>j) by (5ej(>i) = V{ej, ei)SNi (see the small, round, arrows on the e 
axis in figure [3]). In this way, in all the intervals [e^, ej+i], with j ^ i, the A'^- 
and Gj remain unchanged. Therefore, if we maximize the partition function 
with respect to the population of the energy intervals [ej, ej+i], i = 0, 1, . . ., we 
have to take into account the change of the quasiparticle energies due to the 
change of the populations. We shall come back to this method later. 

Another method to calculate the partition function and its maximum, is 
to fix the intervals [eo, ei], . . . , [ej-i, ej], . . . along the e axis. In this case, the 
insertion of the 5Ni particles into the interval [ej,ej+i] changes the values of 
ej (for j > i) given by equation dH) into e^, which should be calculated like in 

e, = e] + Vie'^,e,)Ii + ^ V[e'^, ek]N{ik, h+i) + ^ V[e'j, ei]N{~ek, h+i)- (5) 

A:=0 k=i+l 

Because a{e) is independent of the population, the change of ej(>i) (see the 
5ejS on the left axis in figure [3]) leads to a change in the number of available 
single particle states in the interval [e^-,e^-^^] and therefore a change of the 
number of states in the interval [ej,ej+i] and a change of a{e), according to 
equation ^ . This is the manifestation of FES. 

The exclusion statistics parameters of the FES gas have been calculated 
in Ref. [18j. There we showed in the general case that these parameters obey 
the properties conjectured in Ref. |T5], namely that the mutual exclusion 
statistics parameters are proportional to the dimension of the Hilbert space 
on which they act. Let us take as example a system of bosons of constant 
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Figure 3: The intervals [eo, ei], . . . , [ej-i, ej], . . ., on the e axis, and the cor- 
responding intervals, [eo, ei], . . . , [ej_i, ej], ... on the e axis. The insertion of 
particles in the interval [ei,ej+i] changes the intervals on the e axis, from 
[ci, ej+i], above, if the intervals on the e axis are held fixed. 



interaction potential, V{e, e') = V. Then, the direct and mutual exclusion 
statistics parameters are |18] 

a^i = a^e = Va[e{e)], (6a) 

and 



=9{e-e')V 



d{Ha{e)]} 



de 



a[e{e)]5e = are'6G{e) 



(6b) 



respectively, where 6e is the energy interval in which the mutual statistics is 
manifested and 5G = a[e{e)]6e is the number of energy levels contained in 
it. Having the exclusion statistics parameters ([6]), we can write the integral 
equation (19) of Ref. [T^ for the most probable particle population as: 

[l + n(e)]i-^'^(^) 



/?[/i-e(e)]+ln- 



n(e) 



V / ln[l + n(e')] 



da 
d? 



de'. (7) 



Differentiating both sides of ([7]) with respect to e, we obtain the differential 
equation, 

dn l + Vn{e)a{e) ^de 



de n(e)[l + n(e)] ~^de' 
If we use ^ into ([8]), the differential equation reduces to 

dn{e) 



fi- 



de 



-n(e)[l + n(e)], 



(8) 



(9) 
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which, provided that n(e) is positive and converges to zero at e — > oo, admits 
only the solution 

n{e) = {exp[p{e - ^i')] - (10a) 

This is, of course, the Bose population in the variable e, and relation (jlOap is 
true for any cr(e) and constant V. To determine we plug (jlOap back into 
([7]) and obtain 

H' = IM- VN, (10b) 

where = cr(e)n(e)de is the total number of particles in the system. 
From n(e), the population n[e) follows by a simple change of variable, given 
by equation ([T]). 

Equations (jlOp may seem surprising, but let's have now another perspec- 
tive on the problem. Since F is a constant, the total energy of the system 
reduces to 

E = Y,em + ^, (11) 

i 

which is the energy of the gas in the mean-field approximation. This case is 
easy to treat also from the perspective of the free particle energies, e. If we turn 
back to the division of the energy axis e into the intervals [e^, Ci+i), i = 0, 1, . . ., 
of Ni = //.'"''^ ra(e)cr(e)de particles and Gj = J^^^^ (T(e)de states, we write the 
partition function and we maximize it with respect to the populations n(e), 
then we get 

n[e'(6)]={exp[/3(e'-^)]-l}-\ (12a) 

which is exactly equation (jlOaj) . but with a redefinition of the quasiparticle 
energy, 

e' = e + VN. (12b) 

The quasiparticle energies ()12bp are different from the ones defined in ([T|). 
Moreover, — et is independent of i, unlike ii — ei, which is Z]}=o + 
and depends on i. Another difference is that while the total energy of the 
system can be written as ii^ = hni-, the summation ^[^i gives E+V N"^ /2. 

So, the definition of quasiparticle energies ([I]) leads to the manifestation 
of FES along the quasiparticle energy axis, e. If V is constant, the FES parti- 
cle distribution, ra(e), corresponds to the Bose distribution in the free-particle 
energies, n(e) (equations [TO]) . If a{e) is also constant, then the mutual exclu- 
sion statistics parameters, ajf/ (equation I6bp vanishes and the direct exclusion 
statistics parameters become independent of e: a^i = Va; if Va = 1, then 
the interacting Bose gas may be interpreted as a Fermi gas. From here, the 
thermodynamic equivalence of systems of the same, constant density of states 
and any exclusion statistics [211 E21 EH] follows directly. 
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3 CONCLUSIONS 



This paper is a continuation of Ref. |18j . where we showed that the fractional 
exclusion statistics (FES) is manifested in general in systems of interacting 
particles and the mutual exclusion statistics parameters are proportional to 
the dimension of the Hilbert space on which they manifest. Here we calculated 
the density of quasiparticle states in general and we analyzed in particular a 
system in which the inter-particle interaction potential does not depend on 
the particles quantum numbers (Vij = V). Using the formalism presented in 
Refs. [121 m], we showed that in a gas of interacting bosons, the resulting 
FES quasiparticle population, obtained by the maximization of the partition 
function, is actually the original Bose distribution written as a function of 
the quasiparticle energy, e, instead of the free-particle energy, e. Therefore 
the FES reduces to a change of variable, from e-the free-particle energy-to e- 
the quasiparticle energy~in the population of the single particle energy levels, 
n. Form this result, the thermodynamic equivalence of systems of the same, 
constant, density of states but any exclusion statistics [211 E2l IS] follows 
immediately. 
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